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For distinguishable particles it is well known that Brownian motion and a Feynman-Kac func- 
tional can be used to calculate the path integral (for imaginary times) for a general class of scalar 
potentials. In order to treat identical particles, we exploit the fact that this method separates the 
problem of the potential, dealt with by the Feynman-Kac functional, from the process which gives 
sample paths of a non-interacting system. For motion in 1 dimension, we emphasize that the per- 
mutation symmetry of the identical particles completely determines the domain of Brownian motion 
„Cj , and the appropriate boundary conditions: absorption for fermions, reflection for bosons. Further 

analysis of the sample paths for motion in 3 dimensions allows us to decompose these paths into a 
MH ' superposition of 1-dimensional sample paths. This reduction expresses the propagator (and conse- 

quently the energy and other thermodynamical quantities) in terms of well-behaved 1-dimensional 
fermion and boson diffusion processes and the Feynman-Kac functional. 

> 
On 

i. introduction 

o 

\Q ' In the present paper we introduce a new process that combines multi-dimensional Brownian motion on domains 
with appropriate boundary conditions to solve the many-body problem of fermions or bosons interacting through a 
general class of scalar potentials. This new process is a superposition of orthogonal and independent fermion and 
boson diffusion processes combined in a precise and prescribed way. In combination with the Feynman-Kac functional 
Q this approach allows to write the propagator of the many-body Schrodinger equation as an expectation over the 
| functional along the sample paths generated by the new process. Our method extends the well known techniques 
for quantum models with distinguishable particles to quantum models with identical particles, in such a way 

that the so-called "sign problem" |f7|-|l0f for fermions is solved, and that for bosons and fermions sample paths over 
configurations jll] [ll||8| generated by the permutation symmetry are avoided. This opens the perspective that with 
the proper algorithms the new process would improve the standard approach |l^-|2^1 used in path integral Monte 
Carlo for fermions as well as for bosons. 

Fermion diffusion and boson diffusion are (multidimensional) Brownian motions on a n-dimensional domain D n . If 



o 



X\,X2i ■ ■ -,x n denote the possible components of the positions of the particles on a line, the domain D n is defined by 
the condition x\ > > ■ ■ ■ > x n [p3|. Brownian motion with reflection at the boundary of D n has been studied in 



the context of traffic flow models |24[p5j and will below be identified with the boson diffusion process, because it leads 
indeed to the propagator of n non-interacting bosons on a line. Brownian motion with absorption at the boundary of 
D n has been identified by the present authors |2^,^] as the process which leads to the propagator of n non-interacting 
fermions on a line. We have argued that a 3-dimensional (3D) extension of this process consisting of fermion diffusion 
in one direction and standard Brownian motion in the two other directions is sufficient to obtain the ground state 
energy of an interacting 3D fermion system. The present analysis strongly s upp orts our previous arguments. It should 
be mentioned that our investigation was originated by Korzeniowski et al. |28| and the subsequent discussions on this 
work US 

Because of the fact that for ID problems the stochastic approach to the many-body problem for identical particles 
could be formulated in terms of a fermion or a boson diffusion process, we concentrated our attention to the reduction 
of the sample paths of a 3D process to 3 one-dimensional processes. For the multi-dimensional Brownian motion of 
distinguishable particles this reduction is trivial. Indeed, on the level of the process each 3D sample path is made 
up by ID Brownian motions of the components of the process vector. For identical particles, such a reduction is not 
obvious. But by realizing that the knowledge of the propagator of a non-interacting many-body system of identical 
particles suffices to study how this reduction to independent ID processes can be performed, the character of the 
many-body diffusion could be revealed. The analysis of this reduction is postponed until section 4 and constitutes the 
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main result of the present paper. Before studying this reduction, we first summarize in section 2 the basic underlying 
concepts, borrowed from the stochastic approach to the quantum theory of distinguishable particles. In section 3 the 
fermion and boson diffusion processes are derived and discussed in some detail. 



II. THE FEYNMAN-KAC FUNCTIONAL FOR DISTINGUISHABLE PARTICLES 

In this section we briefly describe how the Feynman-Kac functional and multidimensional Brownian motion can be 
used to solve the Schrodinger equation or to obtain the partition function of a many-body problem with distinguishable 
particles. The basic equations and the mathematical notations are presented. 



A. The 3D configuration space 



For a free particle moving in a 3Z?-space it is well known that the Schrodinger equation becomes a diffusion equation 
if the real time variable is transformed into an imaginary time variable. This means that the motion of the particle can 
be represented by a process |x (t) ;t > j which can take the realization x (t) . The probability density relating 
the position x it) to a previous position x' (to) is Gaussian and given by 



p (x , t; x 1 , t ) 



2-kTi (t - t ) 



3/2 



exp 



vn {x — x') 
' 2h (t - t ) 



(2.1) 



The process |x (t) ; t > j is related to standard Brownian motion |i? (t) : t > j in 3-D with infinitesimal variance 
a 2 = h/m 



dX(t) 



dB (t) 



(2.2) 



If the particle moves in a potential V (x) , the solution of the Schrodinger equatio n ca n be obtained as an average 
over a Feynman-Kac functional. This average is the expectation over the process (2.2). For an arbitrary function 

/ X (t) of the process, the expectation is defined as 



E 3 {/ [X (s)] } = j dx'p (x , t + s; x", t) f (x") . 



(2.3) 



Thus the expectation is the average of / X over all paths that end in the position x at time t + s, starting at 

all possible positions x! at time t. If the function to be averaged is defined by an integral over time, the process is 
constructed in such a way that 



Ex 



X(s) 



% {g [x (s)] } 



ds 



(2.4) 



From the properties of the expectation and the process it is easy to show that for a large class of potentials V (x) the 
solution of the "Schrodinger equation" (in imaginary time) 



is given by 



h^-J> (x, t) = |- V 2 * (x, t)-V (x) * {2, t) 
at 2m 



exp 



V (x,t) = E s ■> f X(t) 



h Jo 



V[X(s))ds 



(2.5) 



(2.6) 



where / X (t) ensures the initial condition on ^ (x, t) . If one wants to calculate the propagator of (|2.5|), this initial 
condition is 
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lim* (x,t) = S (x - x' ) . 

no v ' 



(2.7) 



and hence replacing / X (i) by a ^-function 2.7, one obtains the propagator: 



K{x 1 t\x')=E s \f X(t) 



exp 



-- / V[X(s))ds 



(2.8) 



which solves the Schrodin ger equation (2.5) and satisfies the initial condition (2.7). 

The Brownian motion (2.2) has independent increments in the x-, y- and z-direction. This means that at the level 
of this process there is no difference between the sample paths of one particle moving in 3 dimensions, and 3 particles 
moving each in 1 dimension. At the level of the potential the difference between three ID potentials and one 3D 
potential is obvious. The 3D character of the potential and hence of the problem can (at least partly) be transferred 

to the process by introducing a local drift vector fl X (t) in the process 



dX (t) = fl X (t) 



dt 



dB (t) 



(2.9) 



and the corre spon ding modification in the Feynman-Kac functional J33| , |34[ . But it is clear that changing the process 
according to (2.9) only relabels the state space. 



B. The 3n-dimensional configuration space 



The generalization of the preceding stochastic description of the evolution in the imaginary time domain of a 
quantum system with 3 degrees of freedom to a system with 3n degrees of freedom is straightforward. A point of the 
configuration space given by a 3n-dimensional vector x represents now a state of the process. For convenience in the 
treatment below, we use the following labeling (for brevity writing down the row vector x T ) 



(xx,yi,zi, x 2 ,y2,z 2 , x n ,y n ,z n ). 



(2.10) 



The process {X ( t) : t > 0} with realization x (t) at time t is a configuration obtained according to a straightforward 
extension of (|2.2|): 



(2.11) 



dX (i) = \ —dB (t) , 



where {B (t) ;t > 0} is the 3n-dimensional Brownian motion. The transition probability density for an increment 
(x — x, 1 ) in a time lapse s for this process is given by 



p{ x,t + s;x',t)=(^) 



3n/2 



mix — x'] T \x — x'] 



exp 



2hs 



(2.12) 



The matrix product in ( [2.12 ) can equally well be written as a dot product 



[x — x'] T [x — x'] — (x — x') ■ (x — x') = (x — x') 2 , 



but the matrix notation will be more advantageous for our purposes. 

For the multi-dimensional process {A (t) ;t > 0} the propagator solution of 



* / n 2 ^ d 2 



2m ^-f dx 2 

i=i 1 



can again be written as a Feynman-Kac functional: 

K(x,t\x')=Ejf[X(t)] exp 



V{x) \*(x,t) 



V (X (s)) ds 



(2.13) 



(2.14) 



(2.15) 
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if for / [X (t)] a (5-function is taken to impose the initial condition \imK (x, t\x') = 8 (x — x'). 



It is clear that a straightforward extension of ( [2.9|) with a drift vector ft (X (t)) can facilitate the actual calculations, 
but the formulation of the problem with or without drift is equivalent. Therefore also in the case of 3n degrees of 
freedom, the process ( 2.11 ) that realizes the sample paths for the Feynman-Kac functional does not distinguish 
between n particles moving in 3D or 3n particles moving in ID. This situation will be different in the case of identical 
particles which will be considered below. 



III. PERMUTATIONS AND PROCESSES 



A striking example of how permutation symmetry can be used to simplify joint probability distributions can be found 
in order statistics p^| . The basic idea is that for most samples the order in which the values are measured is irrelevant, 
which means that any permutation of the observed values should have the same probability. This constraint makes the 
joint probability distribution symmetrically dependent. The dependence can be lifted by reducing the sample space to 
ordered sample points: thus for an observation x% , x 2 , . . . , x n the sample point xru , xi 2 -\ , . . . , Xi n \ is considered where 
the observed values are ordered in such a way that xn\ > £(2) > •• • > x ( n )- 

The idea of lifting the dependence implied by the permutation symmetry through a reduction of the domain has 
been put forward before. The basic ingredient of the Bethe Ansatz |37|] is precisely the restriction of the positions 
on a line in such a way that the positions remain ordered. This observation has also been used to solve the nodal 
plane problem for fermions in one dimension || . In particle physics the spatial-temporal distribution of bosons and 
fermions in beams has been obtained using this reduction scheme applied on point processes ||38|-[4l| . For a diffusion 
process we applied an analogous construction to obtain the propagator of interacting fermions Pq^ The process used 
to realize the sample paths was called the "fermion diffusion process" by the present authors. It is a n-dimensional 
diffusion with the appropriate (absorption) boundary conditions. 

We first illustrate the technique for two particles, and summarize the properties of the fermion diffusion process. 
Subsequently it will be shown how this reduction scheme can also be used to construct the n-dimensional boson 
diffusion process. 



A. Two Particles on a Line 



Let X\ and X2 be the coordinates of the first and the second particle respectively. The configuration space is 
two-dimensional: {x\,x 2 ) & R 2 - If the particles are identical the configuration (x±,x 2 ) and the configuration (x2,xi) 
should indicate the same state. For fermions with parallel spin, the anti-symmetry under the interchange of the two 
particles is taken into account by the propagator 

x'i,%2) = p(xi,t;x{,0) p(x 2 ,t;x' 2 ,0) (3.1) 
-p(xi,t;x 2 ,0)p(x 1 ,t;x' 2 ,0). 

We consider this formula for x\ > x 2 and x[ > x' 2 , thus for the position elements of the state space being (xi, x 2 ) G D 2 
and (x' 1 ,x' 2 ) £ D 2 . The boundary of D 2 is defined by x\ = x 2 and denoted by dD 2 . 

It should be noted that the usual factor 1/2! in front of the propagator defined over the configuration space does 
not occur if the motion is restricted to the domain D 2 . Moreover, the propagator on D 2 has all the properties of a 
transition probability density of a diffusion process with absorbing boundary conditions: (i) it is positive on D 2 , (ii) 
it conserves the probability flux if the boundary state is explicitly introduced E3-0|, and (iii) it has the semigroup 
property. These requirements can be checked by direct calculation. The positivity of the propagator can also be 
understood on the basis of a simple geometrical argument: the distance between (£1,0:2) and (x' 1; x' 2 ) both in D 2 is 
always smaller than the distance between (xi,x 2 ) and (x^x^) except at the boundary dD 2 where both distances are 
equal. For a Brownian motion this relation between the distances implies that the probability density to go from a 
point (xi, x 2 ) in D 2 to another point (x' l7 x' 2 ) in D 2 in a fixed time lapse t is always larger than the probability density 
to go from the same point {x\,x 2 ) to the reflected point (x^x^) outside D 2 in the same time lapse. 

Having found the diffusion process for two fermions as a two-dimensional diffusion process on D 2 with absorption 
on the boundary -the so-called fermion diffusion process- it is easy to show that similar considerations hold for two 
bosons, starting from the propagator 

x[,x' 2 \ = p(x- L ,t;x' 1 ,0)p(x 2 ,t;x' 2 ,0) . (3.2) 
+p (xi, t; x 2 , 0) p (xi, t; x' 2l 0) 



XU X2 



-tH/h 



Xi,X 2 



-tH/h 
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Again the propagator ( ]3.2| ) is a transition probability density of a two-dimensional diffusion process on D 2 , but now 
with reflecting boundary conditions on dD 2 . Similarly as for fermions, the required conditions for such a transition 
probability density are easily verified. A nice consequence of the restriction of the state space to the domain D 2 is 
that (for both fermions and bosons) the propagator in D 2 satisfies the initial condition 



lim ( X\ , x 2 
tj.o 



-tH/h 



x[, x' 2 \ = 5 (x x - x[) 5 (x 2 - x' 2 ) if 



(x 1 ,x 2 ) G D 2 
(x' l7 x' 2 ) £ D 2 



It is clear that (|3.l|) and (3.2) can be written respectively as a determinant and a permanent 



X\,X 2 



-tH/h 



X-\ . Xr 



p(x 1 ,t;x' 1 ,0) p(xi,t;x' 2 ,0) 
p(xi,t;x' 2 ,0) p(x 2 ,t;x' 2 ,0) 



(3.3) 



(3.4) 



where £ = +1 refers to a permanent (bosons) and £ = — 1 means a determinant (fermions). This observation allows to 
generalize the process for two identical particles to a process for n identical particles moving in I dimension, because 
the form of the transition probability density will automatically take the boundary conditions into account. The fact 
that the transition probability density has to be zero at the boundary implies absorption at the boundary for the 
fermion problem For the boson problem the transition probability density has to be an extremum at the 

boundary; this implies that its normal derivative at the boundary is zero and as a consequence the boundary condition 
is a reflection j2^^,^]. (Strictly speaking, the boundary dD 2 does not belong to the state space D 2 for fermions. 
But since the propagator is zero on this boundary, there is no point in making the distinction between the state space 
of fermions and bosons as long as the conditions are properly taken into account.) 



B. The Fermion Diffusion Process 



In the preceding subsection it was found that the state space D 2 for two indistinguishable particles is found by 
imposing an ordering on the configuration space, which introduces an additional boundary. The boundary condition 
determines the boson or fermion character of the particles because reflection at the boundary leads to a symmetric 
propagator in configuration space, whereas absorption implies an antisymmetric propagator in configuration space. 
If one considers n particles, moving freely on a line, one obtains (n — 1) extra boundary conditions according to the 
rule 

xi > x 2 > ■ ■ ■ > x n (xi,x 2 , ■ ■ -,x n ) G D n (3.5) 

Defining a diffusion process on D n with the appropriate boundary condition, we ensure that the Feynman-Kac 
functional can be used to incorporate interactions between the particles 0]. Once this is realized, all attention can 
be given to the process. This means that one has to define the transition probability density to go from an element 
of D n to another element of D n in one time lapse. This transition probability density of course has to satisfy the 
conditions for a diffusion process in order to take advantage of the theory developed for multi-dimensional diffusion 
on a domain p5[ . Knowing the transition probability density, one can construct the sample paths using stationary 
and independent increments. 

For the fermion diffusion process, this transition probability density is a Slater determinant, the elements of which 
are the single-particle propagators jj6|. The transition probability density to go from x' £ D n to x G D n in a time 
interval t is 

Pf (x, t; x' , 0) = det \p (xi, t; x'j, 0) | (3-6) 

where p (xi, t; x'j, 0) is the I -dimensional version of the 3-dimensional propagator ( |2.1| ) 
is a transition probability density for a Markov process on D n . This new process |a (t) 
the n-dimensional process {X (t) ; t > 0} given by fl2.11| ) as follows 

X(t) = (f (*) iovt<r dDn 

y > \A(t 9d J for t>r dDn y ' 

The Markov time Tgn n is the first exit time of the domain D n . It should be emphasized that the algorithm described 
in |2(| to realize the sample paths is based on the present formulation of the process. Another realization of this 
process can be obtained with the rejection technique p7| . 



In [ p6| it was shown that pp 
: t > j can be obtained from 
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The semigroup property of the transition probability density pp follows by using the resolution of the unity operator 
1 = J D \x) (x\ dx, whereas the initial condition 

limp F (x,t;x',0) =8{x-x') (3.8) 
follows from the corresponding property for the one-dimensional propagators. 



C. The Boson Diffusion Process 

Like in the case of fermions on a line, the indistinguishability of bosons leads to a state space D n defined in (|3.5| ). 
This means that (n — 1) boundary conditions have to be given to find a process for bosons analogous to fermions. 
The idea that permutation symmetry can imply an ordering for bosons identical to that for fermions was recently 
put forward p6| as a consequence of the arbitrariness of the connection between the statistics of the particles and 
the algebraic properties of the second-quantization operators. Accepting the state space D n , the analysis of D2 for 
bosons suggests that reflection at the boundary dD n of D n leads to a diffusion process that reflects the Bose-Einstein 
statistics. 

Let x and y be two elements of D n , and construct the following permanent 

Pb (x,t;y,0) = perm \p (xi,t; y d , 0)| . (3.9) 
It is clear that is positive for all (x, y) pairs, and that it also satisfies the required initial condition 

]imp B (x,t;y,0) = S(x-y). (3.10) 

Furthermore, in order that p B can be used as a transition probability density it has to satisfy the conservation of 
probability and the semigroup property 



D„ 



Pb (x, i; y, 0) dy = 1; (3.11) 
p B (x, t; y, 0) p B (y, s; z, 0) dy = p B (x, t + s; z , 0) . (3.12) 



The conservation of probability can be derived using the property that a permanent is invariant under an interchange 
of two rows or columns. Hence 

/ pB(x,t;y,0)dy= I —y2,p B {x,t)y p ,Q)dy = —[ peim\p(xi,t;yj,0)\dy = 1 (3.13) 
Jd„ Jd„ nl p n\ ] Rn 

where use has been made of the fact that p (xi, t] yj, 0) conserves probability. 

The semigroup property follows with an analogous procedure by extending the integration domain D„ to R n using 
the permutation symmetry and subsequently using the semigroup property of the single-particle propagators 

f Dn Pb (x, t; y, s) p B (y, s; z, 0)dy=± f D ^ J2 P Pb {x, t; y p , s) p B {y P , s; z, 0) dy 

= 'h Ir" P erm IPs (xi,t;yj,s)\ x perm \p B (yj,s; z k ,0)\ dy = pcrmlps (x t ,t; z k ,0)\ . 

In the last step, the semigroup property of the one-particle propagators gives rise to n\ identical contributions. 

In order to see how the integration over two permanents leads again to a permanent, the following argument might 
be useful. Denote by \y) a fully symmetrized solution of the Schrodinger equation for free bosons, properly normalized 
jl . The resolution of unity is then given by 



1 = ^/ \y)(V\dy (3.15) 

'Rn 



n 



Denoting by H® the Hamiltonian for the z th free particle, and by H° — Y^7=i ^ nc Hamiltonian for n free non- 
interacting bosons, a diffusion from z 6 D n to x 6 D n is given by 



p B (x,t;z,0) = (x e ~ H "^ R z) = - ■ I (x e - ff "C*-'>/ B y) (y e - H " s ' R z) dy (3.16) 



-H°t/h 



1 



n 



H°(t-s)/h 



-H°s/h 



G 



Reduction of all identical contributions to the preceding integral by permutation symmetry then leads to 



p B (x,t + s; z, 0) 





e -H°t/n 


y) (y 


e -H"s/h 


Z )dy = J D 













Pb 0, t; y, 0) p B (y, s; z, 0) dy 



(3.17) 



Therefore pb (x, t; y, 0) is a transition probability density to go from y to x in a time lapse t for a system of non- 
interacting identical particles with Bose-Einstein statistics. The boundary conditions for this process are determined 
by the behavior of Pb(%, t;y, 0) at the boundary dD n . Because VpB (x, t; y, Q) is zero for x e dD n , ps satisfies 
Neumann boundary conditions, leading to reflection for the process at the boundary (|,^|,|5|,[42] Q . 

In the mathematical literature the relation between the boson diffusion process and Brownian motion is referred to 
as the Skohorod equations |24|]25fl . An account on the relation between Brownian motion and diffusion on a domain 
with reflecting boundary conditions can be found in p3. 



IV. MANY BODY DIFFUSION: A PROCESS FOR IDENTICAL PARTICLES MOVING IN 3D 



In this section we will use the known free-particle density matrix of n fermions and n bosons moving in a 3n 
dimensional configuration space to analyze its reduction to propagators on a state space in such a way that these 
propagators are transition probability densities for the processes discussed in the previous section. The starting point 
of our analysis is the projection of the density matrix of distinguishable particles on a density matrix which has the 
correct symmetry properties under permutation of the particle positions: 

pj (x, t; x', 0) = ^ J2 ?P fe» *; 0) • (4-1) 

p 

The projection operator is a weighted average over all elements p of the permutation group. The weight is the 
character £ p of the representation; i.e. £ p = 1 for bosons, whereas for fermions £ p + = 1 for even permutations p+ and 
= — 1 for odd permutations p_. 
The free-particle density matrix is not a transition probability density: it does not meet the criterion that the 
process is in a single state in the limit t { 0, i.e. limp/ (x, t; x', 0) does not satisfy the required initial condition 

S (x — x') . For fermions there is the additional complication of the sign problem because pi (x, t; x! ', 0) is negative in 
certain regions of the configuration space. These objections against the interpretation of pi (x, t; x' , 0) as a transition 
probability density in the configuration space also apply for motions in one dimension. We removed them in the 
previous section by restricting the motion to the state space D n with the appropriate boundary conditions. 

The questions we have to answer for 3D motion are then: what is the appropriate state space and which boundary 
conditions have to be applied? These questions require a further analysis of pi (x, t; x', 0). 



A. The permutation symmetry 



Consider a vector x in the configuration space R 3n , represented as in (2.10) with the x, y and z component of the 



j particle as the (3j) , (3j + 1) and (3j + 2) component of x. Its permutations x p can be represented as 

x p = [p] x 



(4.2) 



where [p] is a 3n x 3n dimensional matrix with one 3x3 identity matrix on each block row and block column, 



corresponding to each particle. For instance, for 2 particles [p] can take one of both forms 
with / = 



[P] 



i o 
O i 



I 

1 o 



" 1 





o" 




" 





" 





1 





and O = 

















1 














The density matrix for the non-interacting identical particles takes the form 

- » *• 0, = ( 2^) *"* h £ e -» ("S M 1 - if 

p 

which can readily be rewritten as 



(4.3) 
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, _. n . ( m \3"/2 / TO ._ _ _, _ /s \ / 1 +n ( m r i- -A 



(4.4) 



B. Projection on even permutations 



We now separate the even permutations {p+} from the odd permutations {p-} , which can be written as {p-} = 
{rp+}, where r is an element of the permutation group which interchanges two particles. Without loss of generality 
we take the first and the second particle. (The same element r has to be used for all elements {p~}). Using the fact 
that £ p + = 1 for both fermions and bosons, one obtains 



p 



r n _ -/\ { eX P (]T [?+] X ■ X') 

[PlX ' X )-^{H r exp [r] 5 • x') 



(4.5) 



where [r] is a 3n x 3n matrix whose operation is to interchange the coordinates of the first and the second particle. 
Hence [r] only differs from the identity matrix in the block column and the block row corresponding to these particles: 



(4.6) 





ro / o • 


• 01 




i o ■ 


• 


r] = 


O I ■ 


• 




.OOO' 


• / . 
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Note that £ r = — 1 for fermions and £ r = 1 for bosons. Elementary algebra then gives 

b+] x ■ x'^j 

\P+] x ■ x'j 

[p+] x-x r 



5Z^ exp (^ \p]x-x 




(4.7) 



where I denotes the 3n x 3n identity matrix (not to be confused with the 3x3 identity matrix /). Since 

is invariant if any two particles are interchanged, the permutation symmetry properties of the density matrix arc 
determined by those of 



cosh \ £ 
sinhi^ 



I-r 
I-r 



\p + ] x ■ x' for bosons 
[p+] x ■ x! for fermions 



Because 



one readily obtains 



where £i and ^2 are the coordinates of the first and second particle in [p + ] x: 



I-r 




- I 
I 




-1 • 
1 • 
• 


OOO • 






-0 


• 


•• 0- 


r [p 4 


]x- 


x 1 = 


(ft -6) 


• (x[ - 



(4i 



(4.9) 









( x P+ ,i \ 


b+] x = 




■ 


X p+ ,2 











The vector Xj indicates the usual three-dimensional position vector, in contrast to x which is a vector of dimension 
3n, as described above. 
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C. The parity of pi (x,t;x',0) and of its components 



The density matrix of n three-dimensional non-interacting identical particles is given by 

Pl {x,t;x',0) = (^) 3 " /2 e (-#t( £ - £ + £ '-'))^^ e w[^][p+]2-' (4.10) 

p+ 

2 cosh (3 ft - &) • (x[ - 4)) for bosons 
2sinh ^1*1 - £2) ■ (x[ — x 2 )J for fcrmions 

which is the key result of the present analysis which allows to answer the questions above on the state space and its 
boundary conditions. Indeed, the decompositions 

cosha-6 = cosh a x b x cosh a y b y cosh a z b z (4. 11) 

+ cosh a x b x sinh a y b y sinh a z b z 
+ sinh a x b x cosh a y b y sinh a z b z 
+ sinh a x b x sinh a y b y cosh a z b z 



sinh a x b x sirih a y b y sinh a z b z (4-12) 
+ sinh a x b x cosh a y b y cosh a z b z 
+ cosh a x b x sinh a y b y cosh a z b z 
+ cosh a x b x cosh a y b y sinh a z b z 

allow to rewrite the density matrix as a sum of 4 terms pi (x, t; x', 0; £ — 0...3) 

3 

PI (x, t; x , 0) = ^ p 7 (x, t; x', 0; £) (4.13) 

2=0 

in which by convention we associate the summation indices £ as follows with the combinations of given parity with 
respect to the reflection plane orthogonal to the indicated direction 



Parity of pj (x,t;x' ,0;t) for bosons 



index 


£ = 1 = 1 1 = 2 1 = 3 


reflection plane 




_L x 


even even odd odd 




even odd even odd 


1 z 


even odd odd even 



(4.14) 



sinh a ■ b = 



Parity of pj (x, t; x' , 0; £) for fermions 



index 


1=0 1=11=2 1=3 


reflection plane 




_L x 


odd odd even even 




odd even odd even 


1 z 


odd even even odd 



(4.15) 



D. The orthogonality relations 

The difference in symmetry of these contributions pj (x, t; x', 0; £) has important consequences. Consider a function 
/ (x") which is invariant under the permutation of the positions of the particles. It then immediately follows that 

/ dx"pi (x, t; x", 0) / (x") pi (x", t; x', 0) (4.16) 
3 

= £ J dx" Pl (x,t;x",0;£) f (x") PI (x",t;x',0;£') 

e,e'=o 
3 

£=0 
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where use has been made of the different parity of the components of pi (x, t; x', 0) to reduce the double sum into a 
single summation. 

The important consequence is that it is sufficient to analyze each component pi (x, t\ x', 0; £) with a given parity 
individually with respect to the interchange of the particles. It is this property which allows to find a diffusion process 
for each pi (x, t; x' , 0; £) separately, as discussed in the next section. 



E. The state space for many-body diffusion 

We now analyze each component pi (x, t;x',0;£) of the density matrix separately. For a given value of I this 
function, defined on the configuration space, can be obtained from a transition probability density defined on the 
state space = D n ® D n (g) D n because it is a produ ct of the transition probabilities of three independent processes, 
each defined on a D n . In other words, the property ( 4.16| ) for the density matrix of free identical particles allows to 
reduce the configuration space to a much smaller state space with independent fermion or boson processes in each 
direction. 

Let {Xpt (t);t>0} be the process that generates the sample paths for fermions moving in 3D. Then this process 
is given according to the following rule 



X Fe (t) 




(4.17) 



where Xp (t) , Yp (t) and Zp (t) denote fermion diffusion processes in the x, y and z direction. Similarly Xb (t) , 
Yb (t) and Zb (t) are boson diffusion process in the x, y and z direction. For bosons the decomposition is as follows 



X B i (t) = 



£ = 
X B (t) 
Y B {t) 

z B {t) 




(4.18) 



For example, the fermion case with £ = 1 is invariant under even permutations {p+} of the particle coordinates. 
Furthermore, under r (interchange of two particles) it is antisymmetric in the ^-direction and symmetric in the y- 
and z-direction. These symmetry properties allow to restrict the transitions x! — » x to a domain D\ = D n (g> D n Cg> D n 
simultaneously satisfying the conditions 



x e Di 



Xl > x 2 > ■ ■ ■ > x n 

yi > y2 > ■ ■ ■ > y n 

z\ > z 2 > ■ ■ ■ > z n 



(4.19) 



with the boundary condition that pi (x,x';t;£ = 1) is zero if during the transition process the boundary dD n is hit in 
the x-direction, whereas it is symmetric with respect to the boundary dD n in the y- and the z-direction. This means 
that a n-dimensional fermion diffusion process is required in the x-direction, and a n-dimcnsional boson diffusion 
process in both the y- and the z-direction. 



F. The Feynman-Kac functional 



The many-body process {t);t > fjj for fermions or |Xb^ (t) ;t > oj for bosons, defined on the state space 

D^, can be used in the same way as {X (t) ; t > 0} to take the interaction between the identical particles into account 



K, 



K { f.b) (MO = T,l=0 K (F B) e i^t\x') 

(x,t\x>) = E w [f [x {FB)t (r)] exp [-\ /„* V (x (FB)t (s)) ds] } 



(4.20) 



Once the propagator of the interacting system on is found with this prescription, the propagator on the 
configuration space can be obtained from the permutations of the particle indices. The function / X (FB) i (r) 
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expresses again the initial condition. By construction the propagators found on D„ have the correct symmetry with 
respect to interchanging of the particles. From the positivity of the expressions pi (x, t;x',0;£) on D^ it follows that 
K [FB) [x, t\x') is also positive on the same domain. 



V. DISCUSSION AND CONCLUSIONS 



The construction of a state space equipped with a diffusion process to provide sample paths for the Feynman-Kac 
functional of an interacting system of identical particles takes into account their statistics by imposing for each degree 
of freedom the appropriate boundary conditions. The spin- states are left out of the picture by assuming that there 
are no spin-dependent interactions involved and therefore the spin as an additional degree of freedom has not to be 
considered explicitly. Of course the spin degrees of freedom are implicitly present because two identical particles are 
only considered indistinguishable if they are in the same spin state. The combined process for a particle to diffuse 
in the configuration space and to change its spin state would require a combination of diffusion and a point process 
fl7|| . The latter would change the number of indistinguishable particles, and hence the dimension of the state space. 
We have only considered the case with a fixed number of indistinguishable particles in the many body problem. 
Given a fixed number of fermions or bosons that are not allowed to change their spin degrees of freedom, a process 
is constructed that takes into account their statistics. All configurations which only differ from each other by a 
permutation of the coordinates of a particle are represented by the same state. For motion in ID this leads to a the 
fermion diffusion process J26|,^] and a boson diffusion process, on a state space D„ with respectively absorption or 
reflection on the boundary <9D„. 

The motion in 3D combined with the permutation symmetry for n identical objects induces reducible representations 
in the configuration space. Fortunately the quantum behavior of a non-interacting many body system is known and 
could be used to decompose the propagator of such a system into independent boson and fermion processes for each 
orthogonal direction of the motion. This reduction is the main result of our analysis because it allows to formulate 
each process in 3D as 4 combinations of 3 independent II? processes with boundary conditions in the appropriate state 
space. It is the combination of the boundary conditions on D n £g> D„ <g> D n that determines wether the process in 3D is 
suited for fermions or bosons. For example a 31? process for fermions can contain one ID fermion diffusion process and 
two ID boson diffusion processes. If one wants to calculate the ground state energy of interacting fermions the two 
ID boson processes may be replaced by a Brownian motion in 2D because the ground state energy of distinguishable 
particles and the boson ground state energy attain the same value in this limit. This observation explains why we 
found that the process used formerly |2^ , Pq| to study ground state properties of interacting fermions leads to the 
correct ground state energy 27 (except when a special direction is chosen to define the ID state space for fermion 
diffusion which leads to an excited state for symmetry reasons). The part of the 3D process for fermions that embodies 
the three ID fermion diffusion processes will decay faster because the odds to be absorbed are higher. Therefore the 
contribution of this process to ground state properties should be negligible. In the case of bosons one may expect for 
the same reason -absorption implies a faster decay- that the low temperature properties can be deduced from the 3D 
process that reduces to three ID boson processes. 

The main achievement of the present approach for bosons is that the sample paths over the configuration space 
taking into account the permutations |ll],|^] are replaced by sample paths arising from 3 independent ID processes 
with well defined boundary conditions. For fermions the significance of our approach is that it provides the solution 
of the so called sign problem Indeed the variance due to the fact that the propagators are antisymmetric under 
interchange of particle coordinates does not arise in our approach. On the level of functional integration we treat 
bosons and fermions with the same method. This is not the case in formerly proposed approaches based on the 
commutation rules for bosons and anticommutation rules for fermions |^,^9|. Our approach relies fundamentally 
on the application of group theoretical concepts to processes, thereby extending previously developed methods for 
distributions |35| and point processes |4l|j40[ |. 

Summarizing we can say that we introduced the state space together with the process that provides sample paths 
for the Feynman-Kac functional of a many body problem with a fixed number of identical particles in a given spin 
state. The propagator over the configuration space can be obtained by the application of permutations to a linear 
combination of such processes. The boson and fermion diffusion process are fundamental processes in our approach 
and their relation to the standard Brownian motion in a domain is settled by the boundary conditions: absorption 
for the fermion diffusion process and reflection for the boson diffusion process. 
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Erratum Februari 10, 1997: 

In section 4.D. we stated that the double sum in Eq. (4.16) 

/ dx"pi (x, t; x" , 0) / (x") pi (x",t; x' , 0) = ^ / dx" Pl (x, t; x" , 0; I) f {x") Pl (x" , t; x' , 0; £') 

reduces to a single sum if / (x") is invariant with respect to the permutations of the particle positions, where x" is 
a 3n dimensional position vector of n particles in a three-dimensional space. This should be corrected as follows: 
the double sum in Eq. (4.16) reduces to a single sum if / (x") is invariant with respect to the permutations of the 
components of the particle positions. This restricts the class of scalar potentials for which the method presented 
can be applied without modifications, i.e. without taking into account the transitions between pi (x, t; x", 0; £) and 
pi (x", t; x', 0; I') for I ^ £'. Furthermore the identity matrix / should be replaced by —I in the second row of the first 
column of the matrix in Eq. (4.8). 
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